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ABSTRACT 




A phenomenological study is made of the shock formation 
of tektites from a meteorite impacting on the earth. The 
calculation shows that if the meteorite that caused the Ries 
Kessel crater ejected large fragments from the parent meteorite 
at speeds of *30 km/sec, these particles would interact with 
the terrestrial ejecta at the start of the recoil trajectory. 

A hypersonic ablation calculation at 20 km/sec is made for a 
single large fragment by extrapolating Hoshizaki's 
• model for the heat transfer rate at an axisymmetric stagnation 
point. 

The results of these calculations show that this physical 
model can yield the size and distribution of the tektite field 
being considered and also account for many of various 

tektite characteristics. The analysis in this paper is 
limited to only observations, thus all refinements which do 
not clarify the basic physical model have been ignored. 



INTRODUCTION 


The purpose of this paper is to examine a physical model based 
upon a shock wave structure in order to explain the terrestrial 
origin of tektites found in the immediate vicinity of craters. 
Tektites in the Czechoslovanian field (moldavites) are considered 
in particular and with emphasis on their origin at the time of the 
Ries Kessel event. 

A preliminary analysis of the propagation of high intensity 
wave discontinuities in impacting solids is currently underway. 

In order to reduce the problem to its essential physical charac- 
teristics an idealized model was chosen and will be discussed in 
a later communication. In this paper we wish to apply only the 
physical results as hypotheses to ascertain the relevance of that 
model to the facts of meteorites at hand. 

A crystalline solid may be stressed into the plastic zone or 
beyond by the inertial reaction to a severe impact (Short, N.M. , 
1966, J. Geo. Educ., XIV, 4, 149). A shock wave then propagates 
through the material causing severe local supraplastic deformations, 
altering the chemistry along its path until it reaches the boundary 
of the material. At this poiht reflection and a transfer of momen- 
tum occurs which will lead to a high velocity mass ejection from the 
free surface of the solid, provided that this momentum transfer is 
above the rupture strength and response time of the solid. The 
directional energy of the shock wave is then converted to direc- 
tional energy of the ejected fragment, thus little shock wave heat- 
ing occurs to the matter that has been dislocated. 



Sufficient data concerning the mechanics of crater for- 
mation and meteorite disintegration are not available in the 
case of the Ries Kessel, whereas a great deal of information 
has been accrued concerning the Barringer Crater event. 
Therefore * in thie report, data is extrapolated from the Bar- 
ringer Crater to bear on the Ries Kessel with the assumptions 
that the Ries Kessel and the Barringer Crater meteorites are 

i 

both metallic, that the target of the meteorites are similar, 

• and that both meteorites impacted at angles no larger than 
55° from the vertical, as is implied from the relative geo- 
metries of the respective craters. 

The diameter of the Barringer crater is 1300 meters, and 

it§ d§pth i§ 175 meters. Aseerding ts And§ 3 r§ (19§§), the §g§ 

of the meteorite which caused the crater is estimated to be 

g 

20,000 years. Its mass is about 2 x 10 kg and its diameter 
ranged from 60 to 80 meters. Fragments from the Barringer 
Crater meteorite have exhibited characteristics which indicate: 

a) shock induced pressures of 130 kb on the outer surface 
and 750 kb in the interior of the impacting meteorite 
(DeCarli and Jamieson, 1961), 

b) temperatures of 700° K on the outer surface and 1100° K 
in the interior of the impacting meteorite. This temp- 
erature gradient is attributed to a shock wave (Lip- 
schutz and Anders, 1961) , 
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c) Chao, Shoemaker and Madsen (1960) found that the 
sandstone beneath the crater contained coesite, which 
is formed by pressures exceeding 20 kb, 

d) Heymann (1964) relating He-3 content to sample depth 
in the meteorite found that specimens farthest away 
from the crater had been nearest the surface of the 
meteor and were the least shocked, while the specimens 
nearest the crater rim had come from deep within the 
meteorite body and had experienced the greatest shock. 

A missile subjected to an impact sufficiently severe 
enough to stress the material locally to a molten state, induces 
a shock wave which carries the effects of the initial deformation 
throughout the material. The classical Rankine-Hugoniot 
relations are no longer valid because the post shock plasma 
reactions fall into the domain of non-relativistic degenerate 
Fermi-gas (NRF) , in which post shock ionization occurs. 

The condition for statistical degeneracy is that the inter- 
particle spacing be equal to the de Broglie thermal wave length, 
i.e., (Landau and Lifshitz, 1958), 

r = h , (1) 

* y2mkT 



in which T = m e c = 2 x 10 9 °K and is the temperature of a 

2k 

relativistic Dlasma. 
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r = average atomic spacing 

X = Compton wavelength 
-l 1 

= 3.86 x 10 cm 

Since T ■» 1100°K# equation (1) shews r s X Q x lo 3 or 
r £ 3.86 x 10 8 cm for a nonrelativistic degeneracy. 

The average interparticle spacing is given by 



in which Z is the atomic numer (57) and N is the particle 

2 4 3 

density (4.7 x 10 per cm ). 

* Thu 8 

r « — x 10"* 8 s. 3.86 x 10~ 

57 x 4.7 

This proves that the nonrelativistic degeneracy of iron 
will be maintained for interparticle compressions up to a 
factor of *10 3 , or volumetric compressibilities -10 9 . 

The most important reaction to take place for the 
purposes of this analysis is the large post shock compression 
that results with the transfer of energy from the external thermal 
modes to internal modes thus cooling the once molten plasma. The 
analysis of this effect is complicated considerably by mathematics: 
we present the results of only one analysis (Skalafuris 1965) . In 
Figure 1 we plot the post shock compressibility vs. the fraction 
of the shock wave energy that is dissipated into internal modes. 

A poly-electronic crystalline substance such as meteorite 
material which has many degrees of freedom , would provide 
a thermal sink to practically all the shock wave energy# thus 



this internal dissipated ratio should be close to unity. In 
view of this fact, a. post shock compression of a thousand is 
perfectly reasonable in which case the material is far beyond 
the state of classical solid state analysis and well into the 
region of a NRF plasma. Thus we subscribe to Bjorks analysis 
as a realistic microscopic formulation of this problem. 

R. L. Bjork (1961) made a theoretical investigation using 
a "Fermi-Thomas-Dirac" theory in the 10 mb range for the im- 
. pacting material and target. Experimental data was available 
in the lower pressure ranges. Bjork assumed that an iron me- 
teor with a velocity of 30 km/sec and weighing 1.2x 10 7 kg 
impacted normally on the ground (tuff) . Bjork gave the im- 
pacting meteorite the geometry of a right circular cylinder. 
His calculations estimate the crater's size to be 150 meters 
deep and 500 meters in radius. If limestone and sandstone 
were the target instead of tuff, (which is more realistic) , 
the meteor would have to have a mass of about 6.4 x 10 kg. 
Clearly, Bjork' s estimations of the size of the meteor for the 
Barringer Crater are much less than those previously quoted. 

We use Bjork' s data to obtain a lower limit to the temperature, 
pressure, and recoil velocity of the Barringer Crater corrobo- 
rating it with an independent shock wave calculation. It can 
be seen from Table 1 which is derived from Bjork' s data that 
the surface pressures are consistent with Anders' analysis. 



TABLE 1 


t - t Q P m ^nternal) p m ^J ear surface) v^gear surface) 
(msec) (Mb) (Mb) (km/sec) 


0.17 

10 

0 

30 

0.36 

10 

0 

30 

3.44 

0.10 

0.01 

-8 

6.36 

0.20 

0.02 

-6 

9.25 

0.05 

0.01 

-6 

24.60 

0.005 

0 

-5 

61 ,00 

0 

0 

-4 


where t - t^ = time after impact 

P s pressure due to hypersonic shock wave in meteorite 
v = velocity of upper surface of meteorite 
msec s milliseconds 
Mb s megabars 
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Theory 

The diameter of the Ries Kessel is from 28 to 29 km 
and it has a depth of 500 meters (Chao and Shoemaker, 1960) . 

If we assume an impact origin for the Ries Kessel, we can 
scale Bjork's data. Following Murphey and Vortman (1961) 
we assume that the volume of the crater produced is propor- 
tional to the energy of tne meteorite. The Barringer Crater 
and Ries Kessel meteorites are assumed to have impacted with 
nearly the same speed (30 km/sec) . Under this assumption, we 
. then have the following relation: 

i 

ft 

T BC _ m BC _ ^BC ^ 

T RK m RK V RK • 

T ■ kinetic energy of impacting meteorite 

m = mass 

V = volume of crater 
| m BC = 64 x 10 6 kg 

V BC = 0.075 km 3 
= 225 km 3 

Therefore, the estimate drawn from Bjork's data leads 
to a mass of the Ries Kessel meteorite of 2 x 10^ kg. 

To estimate the shock pressure built up in the meteorite, 

which is assumed to have the geometry of a right circular 
cylinder, we use the ratio# 
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in which, 


P x dM 1 /dt 1 1/A ] _ 

P 2 dM 2 /dt 2 1/A 2 

M = linear momentum, 

A = area of impacting surface . 


Since the molecular constituents of the target are the 
same, the relaxation time, t, is such that d^ ~ dt 2 ^ t 
for large momenta interactions. Therefore, 




( 6 ) 


We have already hypothesized that » v 2 (impact velocities) ; 
therefore, let 

W m if f (?) 

in which. 


V = a characteristic function of the meteorite, 
p = viscosity, 

r = geometrical shape factor 

\ 

We define f^ 2 = — Which is independent of radius. 

2 

Theoretical results state that the stress is independent of 
the mass at impacting object, (Clebsch, 1883), while experiments 
tend to imply that it does depend on the mass (Voigt , 1915) . Since 

this is a phenomenological study, we adopt the latter conclusion. 
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# * 


Therefore, 



(8) 


Because the meteorites are assumed similar except for their size, 
f^ 2 » 1; hence, 



(9) 


From Eq. (4) we find that R s radius of meteorite, and 
R 1 

jTm 15 . 

R 2 

To determine the recoil velocity (v^ of the meteor 
fragments, we., use the ratio, 



V v i 



( 10 ) 


in which e = specific energy, and, 



(ID 
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Thus , 



3.9 


( 12 ) 


The temperatures for the Barringer Crater are estimated 
of 1100 °K internal and 700 °K external. Because of the short 
duration of the compression, the impact will be adiabatic with 
the temperature and pressure given by the following law. 




y/y-1 


(13) 


If only translational modes are thermalized, y = 5/3; thus, — 

t 2 

The temperature in the inner region of the Riea Keesel 

meteorite is T ^ 3300 °K and the temperature in the outer 

RKt 

region is T 2100 °K. 

KKq 

Because of the gross energies involved at impact, we assumed 
that the major portion of the meteorite is vaporized and the 
equations of an ideal gas hold. The front outer temperature 
calculated is sufficient to throw off the target rock, at the 
interface between the target and the meteorite, in a molten 
state. 


~3 


In Table 2 we list the pressures, recoil velocity, and the 

time after the impact for the Rios Kessel obtained by scaling 

Bjork's data with our formulas. The data from Table 2 is 
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TABLE 2 


t - t 0 
(msec) 

P (interior) 
(Mb) 

P (exterior) 
(Mb) 

v (rear eurfaee 
(km/sec) 

0.17 

150.0 

0 

30.0 

0.36 

150.0 

0 

30.0 

3.44 

1.5 

0.15 

-31.2 

6.36 

3.0 

0.30 

-23.4 

9.25 

0.75 

0.15 

-23.4 

24.80 

0.075 

0 

-19.5 

61.00 

0 

0 

-15.6 


plotted in Figure 2. From Figure 2 it can be seen that the 

internal and external pressure undergo a single maximum, 

differing by an order of magnitude. The graph also shows that 

the recoil velocity peak occurs behind the internal pressure 

peak and falls off as the external pressure rises. 

From previous calculations, Eq. (9) , we estimated the 

size of the Ries Kessel meteorite to be roughly 300 meters in 

radius. The depth of the fragment blown off of the impacting 

meteorite corresponds to the thickness of the shock wave 
* 

induced on contact. A shock wave converts organized motion to 
random thermal motion; hence the thickness (d) of ejected 
matter can be estimated by the equality (see appendix A) , 

d = zmhv^ , (14) 

f kT 

which, 

z s atomic mass number, 

f s degrees of freedom per heavy particle, 

L s length of the meteorite (300 m) , 
v s impacting velocity (30 km/sec), 
kT h thermal energy (at 0.33 x 10 4 °K) , 

If we choose iron with z = 57 and f = 3 consistent with 
assumption (13) , then we obtain a thickness of about 20 meters 
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We assume that this fragment has the geometry of a hemi-ellip- 

soid in which case its radius would be 90 meters. The volume 

5 3 

is then about 3 x 10 m . 

We assume that when the meteorite impacted, the terrestrial 
material was ejected at the time of the recoil velocity peak 
which was at 3.44 msec. From an examination of data we conclude 

that the fragment blew off after 24 msec (where p ex t e rior = ^ there- 
by intercepting terrestrial material from ground zero to about 400 

meters in the atmosphere. 

It is known that under hypersonic flow conditions the 
projectile is followed by a recirculation region which is 
approximately twice the length of the projectile (Fig. 3) 
thus, the ejecta will fill up the recirculation region which 
is at a relatively low pressure. 

If the meteorite fragment blew off at 24 msec with 
a velocity of *-20 km/sec (scaled from Bjork's data), it will 
intercept this ejected material and carry it on the fragment's 
surface as it goes into its trajectory. 

The question of whether or not a meteorite fragment can 
survive such a severe recoil velocity and still remain intact 
seems to be resolving in favor of a coherent mass, because the 
ejection represents a transfer of the shock wave group velocity 
into direction motion rather than post shock thermal dissipation. 
However, these results are yet premature and we still consider 
this an open question for further research (Remo and Skalafuris, 
1967) . 
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The volume of the ejecta that the fragment intercepts 
can easily be calculated. If 5 m of ejecta lie on the top 
surface of this fragment, (an arbitrary but not unreasonable 
figure) then the ejeated volume is . 

We now look into the aerodynamic effects on this pro- 
jectile, being primarily concerned with the heat flux to this 
blunt body for ablation and the energy expended in traveling 
through the atmosphere (Fig. 4). In calculating the heat 
flux we start with a laminar stagnation point calculation. 

The Knudsen number in this problem remains less than unity 
(due to the large dimension of the body) . Therefore we are 
assured of operating in the region of hydrodynamic flow. 

Hoshizaki (1962) presented a relation for calculating 
the heat transfer rate at an axisymmetric stagnation point 
in the velocity range from 1.8 to 15 Jan per second. The 
equation has the following form: 



We will use this relation (15) in the range of 20km/sec. 
This expression includes the effect of variable Prandtl 

^ C p . D 12 

^ and Lewis numbers and also takes into account 

P 

ionization, from equilibrium chemistry. 


(15) 
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C s specific heat 
P 

k = thermal conductivity 
D^2 h chemical diffusion 
Expression (15) can be simplified (Appendix b) : 



C 

air s specific heat of air 

= radius to center of curvature of projectile 
= velocity of projectile 

t 

= density of medium in front of shock wave 
p = density of medium in stagnation zone 

= pressure of medium in front of shock wave 
P = pressure of medium in stagnation zone 

W 

q = heat transfer rate at the stagnation point 
s 

M s mach number 

oo 


Table 3 represents the amount of heat transferred in the 
stagnation zone, as a function of height over the earth. Ef- 
fects of the heat transferred from ionic interactions can be 
neglected. 

It is apparent from Table 3 that the thermal heating 

takes place mainly in the first 25 km of the trajectory. If 

the fragment comes off at an angle of 37 ° , the total heat 

19 

transfer to the body at a height of 300 km is 5.8 x 10 ergs. 
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TABLE 3 

HEAT TRANSFER AT AERODYNAMIC STAGNATION POINT 
(Densities from Rocket Panel, PHYS . REV. 88, 1027 (1952)) 


height 

Pf 

^s 

Q 

(km) 

(gm/cm 3 ) 

(ergs/cm -sec) 

(ergs/sec) 


0-25 

300 

X 

io ' 6 

5.8 x 

o 

H 

O 

3.5 

X 

io 19 

25 - 50 

7 

X 

10“ 6 

1.3 x 

io 9 . 

7.9 

X 

io 17 

50 - 75 

2.6 

X 

10 “ 7 

3.4 x 

10 8 

2.1 

X 

10 17 

75 -100 

3.9 

X 

10- 9 

— 





100 -125 

1 

X 

10- 10 

— — 





125 -150 

1.18 

X 

10- 11 







where Q = q g A 

A = surface area in the stagnation zone 
A « 6.1 x 10 4 m 2 


Heat transfer in stagnation zone 
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From this analysis the loss in velocity due to the heat tran- 
sfer in the stagnation region is- 0.1 km/sec. If the heat of 
fusion-vaporization is about 10^ ergs/gm, we can melt off 

about 2 x 10 gms of ejecta. At 300 km above the earth's 
surface, the loss of velocity due to the gravitational potential 
is 2 km/sec. 

We now consider the loss of energy due to aerodynamic 
drag. Consider the pressures in the compression zone and the 
stagnation zone. 

From Appendix Eq. (3) , 


P s = P f + 1/2p f U f 2 ' 


and 


(17) 


in which, 



( 18 ) 


= recirculation zone pressure, 
P Q = drag pressure , 

= 328 x 10~^ gm/cm 2 , 

Uf “ 25 x 10 7 cm/sec , 
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M 2 molecular weight = 29 gm/ruole , 

R = universal gas constant = 

8.315 x 10 7 erg-deg*"^ c-mole*"^ , 
T = absolute temperature « 230° K , 

= 22 x 10 4 ergs/cm 2 , 

P s ^l/2 Pf U f 2 , 

(19) 

9 2 

P m 10 gm/cm - sec . 
s 


'We then have, 


F = / P • dA , 

o 


F 2.4 x 10 17 gm/cm - sec^ f 


( 20 ) 


in which, 

dA i effective front normal area , 

F = force acting on the compression front of 

the projectile . 

• ' 11 

Now • 


P b - l/K (P s > . 

in which K = constant fraction. Then , 


(21) 


P 


drag 



(1 - l/K) . 


( 22 ) 
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Recent observations in satellite entry and shock tunnel 
testing have shown that while the pressure in the recircula- 
tion region is much less than the pressure in the compressional 
region, the recirculation region pressure is not negligible. 

It is also apparent that the majority of the drag occurs in 
the first 50 Ian of flight, because falls off very rapidly. 
Taking into consideration back pressure (P^) in the recircula- 
tion region, we can estimate p drag by taking K * 10 and using 
Eq. (22) . The total drag force acting on the projectile which 
comes off at a 37° angle causes a velocity change of about 
5 km/sec on the fragment. At a height of 300 km the speed of 
the fragment is reduced to about 13 km/sec. Escape velocity 
at this height is about 9.2 km/sec. 

Though it is clear from the previous results that the 
fragment will escape from the earth, some material on the 
surface of the meteorite will also return, because the molten 
ejecta runs down the side of the fragment and approaches its 
base, where it undergoes a severe pressure gradient caused by 
expansion waves emanating from the fragment base as shown in 
Fig. 5. Since the viscous force r is small, this pressure 
gradient will accelerate this virtually unbound stream of 
molten drops from velocities of zero to almost the free stream 


velocity U^. This glass-like stream will move in a direction 
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opposite to that of the fragment, and will have a velocity "v" 
relative to the earth such that 0 < v < 13 km/sec. Therefore, 
some droplets escape from the earth and others fall back 
covering a long narrow area at a distance from the Ries Kessel 
of roughly between 50 and 2000 km. Passfe 3EI &£ this analysis 
will investigate in greater detail the dynamics of the stream 
of confined droplets as they pass over the meteorite surface, 
and leave the meteorite to enter the atmosphere. 

The dense spray of molten droplets will not be dispersed 
upon leaving the meteorite fragment because the atmosphere 
above 25 km is very thin and can cause only negligible drag 
effects, and so droplets can even become rigid before entering 
the atmosphere where they will undergo aerodynamic ablation 
(Chapman, et.al., 1962). Since the fragment giving off these 
droplets is considered to be moving in a straight line, acting 
as- a line-source, thereby accounting for the distribution of 
the moldavite field. 


Conclusion 

This model accounts for the observed assymetry in the 
tektite field about the Ries Kessel. It also accounts for its 
well defined boundaries since there is only a limited velocity 
range with which tektites could leave the meteorite fragment. 
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The dimension of the Lake Bosumtwi Crater on the Ivory 
Coast which is 10.5 kilometers in diamter (Vand, 1965), yields 
a mass of 3.4 x 10^ kg, large enough to fit into the above 
model. 

It ie interesting to note that the calculated maeooe of the 
bodies that have impacted at the Ries Kessel and Lake Bosumtwi 
sites are of the order of lO^kg. These large masses suggest 
that the impacting projectiles are asteroids. The dates of the 
.impact of these projectiles and their implications have been 
discussed (Cohen, 1963; Fleischer and Price, 1965; Faul, 1966; 
Zahringer, 1963; Barnes, 1963; Gentner and Zahringer, 1959; 
Schnetzler et.al, 1966; Gentner et al, 1961, 1963; O'Keefe, 1963). 

We stress that in this calculation, only numbers were 
derived which were essential to understanding the basic physical 
model. Many parameters such as recoil velocity, meteorite size, 
fragment size, could be changed without sacrifice to the results 
of the basic model. 
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APPENDIX A 

The kinetic energy of organized motion per particle in the 
meteorite before impact is equal to, 

E==»3ALzmv^, (1) 

o 

in which, 

A = cross sectional area impacting / 

L = length of the object , 
z = atomic mass / 
m E mass of a proton, 
v = velocity . 

This energy is delivered to translational modes whose 
collective capacity is , 

E t = -2 (A d) kT , (2) 

in which, 

f = degrees of freedom of the atomic species , 

d = thickness of the region of translational 
relaxation of the shock wave, 

k z Boltzmann constant, 

T 5 temperature . 



Equating the two quantities yields, 


d 


z 


= 7 


_ 2 
mLv 

“kT"" 


where d is the thickness of the ejected matter. 
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APPENDIX B 


The stagnation point heat rate is proportional to 
i / 2 dU e 

, where is the velocity gradient. An 

dU 

expression for e is 

ds 



dU* 

ds 


1 _ 




( 1 ) 


Conservation of momentum for a normal shock is 
represented by the following equation: 


P f + 1/2 p f U f 2 = P s + 1/2 p s U 3 2 (2) 

There exists a low mach number behind the shock 
wave in the stagnation zone. Therefore, 



1/2 Pf U f 2 


(3) 


Then for the velocity gradient we have 



dU^ 

ds 


( 4 ) 
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Where 1/6 > if > 1/10. 
p s 

Since the heat in front of the stagnation region is much 
greater than the heat on the surface of the projectile. 



Then 



axr 



1/4 


M 3 - 19 

10 4 j 



(5) 


Hoshizaki computed the heat transfer for an equili- 
brium boundary layer by using Hansen’s equilibrium 

transport properties. For flight velocities between 

N 

1.8 and 15 km/sec the heat transfer parameter ja has 

K 

the following temperature and flight velocity dependence 



0.2 


0.478 


w 


U, 


-0.31 


500° K 


10 


(6) 


where T w = temperature at projectile wall 

N u = Nusselt number 
R e = Reynolds number 



25 km/sec 


0.235 
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The stagnation zone extends approximately 30° on 
either side of the horizontal (Fig. 3) . 

tfhe amount of heat transfer is calculated, with 
fciw sallowing values tar the paaeamaters i 


= 13.2 m 

C air = 7.2 x XO 12 Z™ 1/2 , 

gin - sec^ 

= 25 km/sec 
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FIGURE CAPTIONS 


Fig. 1 - Post Shock Dissipational Compressibility 

Fig. 2 - Interior, exterior pressures and recoil 
velocity plotted as a function of time 
after impact. 

Fig. 3 - Shock effect due to a projectile traveling 

at a hypersonic speed through the atmosphere. 

Fig. 4 - Diagram showing the heat transfer rate at 
an axi symmetric stagnation point. 


Fig. 5 - Effect of the pressure gradient at the base 
of the projectile on the ablated material. 
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